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abstract 


Progress  on  both  the  short- and  long-period  constituents  of  the  International 
Seismic  Month  is  reported.  Further  work  on  understanding  NOUSAR's 
complexities  is  described.  Evidence  is  given  of  the  impact  of  lateral  het¬ 
erogeneity  within  the  earth  on  seismic  observations.  We  describe  our  first 
work  on  evasion  and  counter-evasion  problems. 
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SI  MM  ARY 


This  is  the  nineteenth  Semiannual  Technical  Summary  report  of  M.l.T.  Lincoln  Laboratory's 
Seismic  Discrimination  Ciroup.  Our  work  is  directed  toward  the  effective  use  of  all  available 
seismological  data  in  determining  the  Lnited  States'  capability  to  detect  and  identify  underground 
nuclear  weapons  tests.  For  several  years  a  technique  has  existed  which  successfully  identifies 
a  large  majority  of  seismic  events  of  magnitude  greater  than  a  certain  value.  Our  concern  at 
present  is  to  i inti  ways  to  lower  this  threshold,  to  understand  occasions  when  the  discriminant 
fails  to  work,  and  to  identify  features  within  the  earth  which  can  degrade  our  detection  and  dis¬ 
crimination  capability. 

Lor  the  past  18  months,  we  have  been  concentrating  attention  on  the  capability  of  deployed 
seismic  instruments  during  a  period  in  February  -  March  1972  known  as  the  International 
Seismic  Month  (ISM).  Wc  started  with  an  attempt  to  produce  a  comprehensive  epicenter  list 
based  on  short-period  detections  at  l.ASA,  MORS  A  It,  selected  World-Wide  Net  stations,  Cana¬ 
dian  Network  stations,  and  Swedish  and  f.  K. -sponsored  arrays.  This  work  required  the  devel¬ 
opment  of  a  major  new  set  of  computer  programs  called  DADS  (Data  Analysis  and  Display 
System).  With  the  aid  of  this  interactive  facility,  the  epicenter  list  is  near  completion;  a  poi  - 
tion  ol  the  liual  list  is  shown.  In  order  to  produce  a  satisfactory  location  which  takes  into 
account  not  only  arrival  times  of  compressional  phases  but  also  any  dT/dA  information  that 
arrays  may  supply,  a  completely  new  algorithm  has  been  developed  and  has  been  successfully 
used  in  reducing  data  from  the  ISM.  More  recently,  we  have  started  to  collate  long-period 
information  which  comes  from  arrays,  conventional  single  stations,  and  the  l.amont-sponsored 
long-period  experiment.  Results  from  this  side  of  the  month's  study  will  follow  rapidly  on 
completion  of  the  epicenter  list. 

The  Norwegian  array,  NORSAll,  continues  to  be  the  subject  of  much  of  our  attention.  The 
puzzling  complexity  of  short -period  P-wave  observations  across  the  array  is  further  examined, 
and  order  appears  to  be  emerging  from  some  of  the  complicated  recordings  of  Eastern  Kazakh 
events.  First  results  of  an  "array  diagram"  approach  to  NORSAR's  mislocation  anomalies  are 
given  and  they  show  features  very  similar  to  those  observed  at  l.ASA.  Multipath  observations 
of  surface  waves  are  described  -  the  location  of  NO  USA  R  may  make  it  less  vulnerable  to 
multipathing. 

We  continue  to  have  a  strong  interest  in  the  effects  of  heterogeneity  in  the  earth  on  seismic 
signals.  A  model  of  the  crust  and  upper  mantle  in  the  western  l  nited  States  has  been  used  for 
ray  tracing  of  surface  waves  in  order  to  give  some  idea  of  the  elTect  that  heterogeneity  within 
a  continental  region  has  on  l.ASA 's  surface  wave  reporting  ability.  Further  evidence  for  lateral 
hetei  ogeneily  deep  within  the  earth  is  presented.  Multipathed  body  waves  are  reported  from  an 
event  distant  95  from  l.ASA,  and  an  amplitude  anomaly  is  described  for  events  in  the  Pacific 
recorded  at  l.ASA.  In  both  cases,  some  sort  of  heterogeneity  is  necessary  in  the  middle  or 
deep  mantle  to  explain  the  observed  phenomena.  An  analysis  oT  l.ASA's  array  diagram,  in 
terms  of  potential  constituents  from  different  depths,  is  given. 

We  are  devoting  an  increased  amount  of  our  effort  to  problems  of  evasion  and  counter¬ 
evasion.  Two  specific  evasive  techniques  are  considered:  multiple  shots,  and  hide-in-earthquake. 
Seismograms  have  been  synthesized  to  give  an  idea  of  what  would  be  needed  to  mimic  complexity 
and  polarity  reversals  of  the  depth  phase.  A  seismogram  which  distorts  the  M^m^  relationship 
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SEISMIC  DISCRIMINATION 


I.  INTERNATIONAL  SEISMIC  MONTH 

A.  EPICENTER  LIST  PREPARATION 

International  Seismic  Month  (ISM)  activity  has  been  focused  on  two  areas  during  the  past 
several  months.  First,  it  was  decided  to  review  all  short-period  data  and  issue  a  revised  and 
more-documented  event  list.  Secondly,  the  collection  and  organization  of  long-period  observa¬ 
tions  were  started  in  earnest.  By  30  June  1973,  new  event  lists  and  long-period  observations 
will  have  been  sent  to  all  active  participants. 

Previous  issues  of  our  event  list  were  deficient  in  several  respects.  Hypocenter  informa¬ 
tion  was  from  several  unspecified  and  heterogeneous  sources.  No  quantitative  estimates  of  ep¬ 
icenter  quality  were  supplied.  Similar  comments  could  be  made  about  the  body  wave  magnitudes 
supplied.  The  new  list  is  considerably  improved  in  these  respects. 

To  facilitate  generation  and  future  manipulation  of  event  lists,  several  major  new  programs 
were  written  to  operate  in  the  Data  Analysis  and  Display  System  (DADS)  environment.  Notable 
among  these  are  /RES  and  /LOC.  These  two  are  programs  to  calculate  residuals  and  to  locate 
events.  The  core  of  these  programs  is  advanced  subprograms  for  residual  calculation  and  loca¬ 
tion.  The  locate  subprogram  uses  azimuth  and  /elocity  observations  as  well  as  arrival  times. 

In  addition,  it  can  make  use  of  most  of  the  major  seismic  phases  so  that  no  information  is  dis¬ 
carded.  Finally,  it  supplies  a  complete  four-dimensional  space-time  error  analysis  from  which 
hypocenter  and  epicenter  quality  can  be  obtained.  More  details  concerning  the  subprograms 
follow. 

/RES  and  /I..OC  work  directly  with  large  data  matrices,  with  each  row  representing  an  ar¬ 
rival.  Arrivals  from  many  events  are  mixed  together.  The  data  matrix  contains  a  column, 
GRPNUM,  used  to  indicate  which  arrivals  all  go  with  a  single  event  by  assigning  the  same  value 
to  GRPNUM  for  all  of  the  arrivals  for  a  single  event.  /RES  will  calculate  time  and  velocity  re¬ 
siduals  and  enter  them  into  the  data,  along  with  the  hypocenter  used.  As  arrivals  are  added  or 
removed,  phase  assignments  changed,  event  associations  changed,  etc.,  it  is  the  task  of  the 
DADS  data  structure  and  /RES  to  be  certain  that  all  data  remain  internally  sell-consistent.  Hy- 
pocenters  may  be  input  by  cards  or  obtained  internally  from  the  DADS  data  base  itself. 

/l.OC  is  similar  to  /RES  in  many  ways,  but  actually  locates  the  event  as  well.  The  hypo- 
center  may  be  automatically  accepted  and  many  events  processed  in  batch,  or  the  user  can  ob¬ 
tain  a  display  preliminary  to  accepting  the  location.  The  display  shows  various  quality  factors 
including  the  residuals  for  the  starting  and  final  hypocenter.  The  user  may  modify  parameters, 
remove  arrivals,  restrain  depth,  etc.,  and  rerun  the  event.  Alternatively,  he  may  decide  to 
reject  it  totally  or  accept  it. 

During  the  past  several  months,  all  ISM  events  have  been  processed  by  /RES.  The  results 
were  reviewed  and  corrections  made  to  the  data.  Each  event  was  then  processed  by  /LOC  using 
the  corrected  data.  Tno  results  were  again  reviewed,  and  /l.OC  was  rerun  on  any  events  having 
further  changes.  All  events  on  our  new  list  will  have  been  located  by  /LOC.  Events  still  causing 
difficulty  or  with  insufficient  data  will  not  be  included. 
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The  new  event  list  will  include,  in  addition  to  hypocenter  and  origin  time: 

(1)  The  largest  semiaxes  of  the  epicenter  error  ellipse  and  of  the  hypocenter  error 
ellipsoid. 

(2)  The  number  ot  arrival  times  and  velocities  used  in  location,  and  the  number  asso¬ 
ciated  with  the  event  but  not  used. 

(3'  An  value  calculated  from  amplitude  and  period  data  contained  in  DADS,  exelud¬ 
ing  reports  more  than  0.5  unit  away  from  the  given  value. 

(T  The  number  of  values  averaged  to  calculate  mb  and  the  number  excluded  by  virtue 
of  being  too  far  from  the  average 

(5)  Minimum  and  maximum  station  m^  values. 

Detailed  corroborating  information  on  any  event  will  be  available  to  participants  on  request. 

Table  1-1  is  a  short  sample  of  the  kind  of  event  list  which  will  be  distributed. 

I  p  to  the  present,  ISM  activity  with  short-period  data  has  been  limited  to  data  organization 
and  event  list  preparation,  as  outlined  above  and  in  previous  SATS.  However,  Fig.I-1  shows 
our  epicenter  quality  as  a  function  of  our  mb  values  for  events  during  a  three-day  period.  The 
quality  factor  (MAXAX2)  shown  is  the  maximum  principal  semiaxis  (km)  of  the  epicenter  error 
ellipse.  The  exact  interpretation  will  require  further  research,  but  it  is  probably  very  conserva¬ 
tive  to  assert  that  the  true  epicenter  will  be  no  more  than  MAXAX2  away  from  the  calculated  ep¬ 
icenter  with  high  probability.  Note  that  quality  is  excellent  down  to  about  mb  =  4.5,  and  then 
deteriorates  quite  rapidly.  It  should  be  noted  that  the  event  at  mb  =  4.8  with  poor  quality  actually 
has  mb  determined  from  only  one  measurement  and  that  there  is  poor  distance  and  azimuth  con¬ 
trol.  An  important  point  to  note  is  that  one  can  lie  automatically  informed  about  potential  errors 
in  calculated  epicenters.  Clearly,  interpretation  of  the  ISM  data  for  detection  and  location  will 
be  a  considerable  challenge. 

Long-period  readings  and  measurements  have  been  accumulated  from  several  sources  and 
organized  into  an  integrated  data  set.  Final  event  associations  and  analysis  have  not  been  under¬ 
taken,  although  preliminary  associations  by  outside  participants  and  Lincoln  are  indicated  in  the 
data  base.  Several  thousand  long-period  vertical  measurements  of  signal  and  noise  were  con¬ 
tributed  by  Lincoln  from  visual  analysis  of  some  15  WWSSN  stations  for  the  entire  time  period. 
The  Seismic  Array  Analysis  Center  did  frequency-wavenumber  analysis  of  selected  time  periods 
for  the  vertical  component  of  the  large  arrays,  and  supplied  the  computer  output  to  Lincoln. 

These  data  have  been  sifted  by  us  and  entered  into  the  data  base  for  further  interpretation.  Inso¬ 
far  as  possible,  basic  observational  data  from  all  other  participants  have  been  included  as  fur¬ 
nished  by  them,  but  derived  quantities  (such  as  Mg  values)  have  not. 

R.  T.  Lacoss 
R. M.  Sheppard 
R.  E.  Needham 

B.  DATA  ANALYSIS  AND  DISPLAY  SYSTEM  (DADS) 

The  DADS  continues  to  be  refined  and  expanded.  Some  significant  features  that  have  been 
added  are  the  /SEARC  and  /PLOT  programs  which  are  described  below,  and  programs  to  com¬ 
pute  source  locations  and  residuals  for  the  ISM  which  are  described  in  Sec.  A  above. 

DADS  is  an  interactive  display  system  designed  to  operate  on  tabular  data.  Generally,  the 
columns  are  referred  to  as  parameters  and  have  names  assigned  to  them.  The  rows  are  events 


or  arrivals  and  are  numbered  sequentially.  A  selector  is  a  named  column  that  has  the  special 
function  of  flagging  certain  of  the  events  that  form  a  meaningful  subset  of  the  data. 

Previously,  selectors  could  only  be  created  by  specifying  the  events  by  number  or  by  writing 
a  special-purpose  Fortran  program  to  choose  them.  Now,  the  user  merely  types  in  a  list  of  con¬ 
ditions;  and  /SEA  111.’  sets  up  a  selector  containing  all  the  events  which  meet  those  conditions.  For 
example,  the  user  might  want  to  set  up  a  selector  B1(!KV  which  consists  of  events  where  the  mag¬ 
nitude  parameter  (Alii)  lias  a  value  greater  than  5.0  and  the  station  reporting  the  arrival  (A.STA) 
is  either  1ASA  or  NOllSAR,  The  command  would  lie; 


/SKA  11C 

1.  A.STA  KQ  I  AO  NAO 

2.  MB  CiT  5,0 
SKI  ,=BlCi  KV  1  2 


The  program  stores  the  conditions  typed  in  and  displays  them  on  the  scope.  The  user  can 
thus  set  up  several  selectors  with  different  combinations  of  the  same  basic  conditions  without 
having  to  retype  the  conditions  each  time.  A  search  imposing  several  conditions  on  up  to  8000 
events  executes  in  less  than  1  sec. 

The  /PLOT  program  produces  a  scatter  plot  of  any  two  parameters  for  the  events  in  a  given 
selector.  Scaling  may  be  automatic  or  preset  by  the  user.  The  user  may  optionally  control  the 
character  plotted  for  each  event  by  previously  setting  up  a  parameter  column  which  contains  that 
character  for  each  event.  Figure  1-1  is  an  example  of  the  output  of  this  program. 


L.  J.  Tarek 


C.  KXTKNSION  Oh’  STANDARD  13V ENT  LOCATION  PROCKDUllES 


Conventional  methods  for  calculating  seismic  hypoccnter  locations  make  use  of  only  the  ob¬ 
served  arrival  times  of  the  P  and  PKP  phases.  This  is  not  too  serious  a  restriction  for  large 
events  which  are  widely  recorded,  hut  for  small  events  it  frequently  happens  that  potentially  val¬ 
uable  data  are  discarded.  Array  measurements  of  dt/dA  and  azimuth,  for  example,  cannot  be 
included.  Nor  can  arrival  times  of  other  phases  (c.g,,  S,  PKKP,  etc.)  which  often  contain  signif¬ 
icant  information,  especially  on  focal  depth  and  origin  time.  The  standard  iterative  least-squares 
location  method,  however,  is  easily  adaptable  to  include  such  data. 

For  use  in  processing  data  for  the  ISM,  we  have  developed  a  generalized  location  program 
which  accepts  any  combination  of  dt/dA-azimuth  and/or  arrival  time  data  for  any  phases  (for 
which  travel  time  tables  are  available).  At  present,  the  phases  in  use  are,  P,  PP,  PcP,  PKP 
(A-B,  B-C,  and  D-F  branches),  PKKP  (A-B,  B-C,  and  D-F),  P'P'  (A-B,  B-C,  and  D-F),  S, 

SS,  ScS,  SKS  (A-C,  C-D,  and  D-F),  including  depth  phases  in  all  cases. 


Another  deficiency  of  standard  location  methods  lies  in  the  technique  used  to  estimate  the 
accuracy  of  calculated  hypoccnters.  The  magnitude  of  the  errors  of  observation  is  estimated 
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In  fact,  however,  these  "errors"  are  produced  primarily  by  lateral  velocity  variations  in  the 
earth  and  are  thus  systematically  related  to  each  other  and  can  introduce  systematic  errors  into 
the  calculated  locations.  In  other  words,  velocity  variations  in  the  earth  produce  travel  time 
anomalies  which  often  can  be  "canceled  out"  by  a  suitable  mislocation  of  the  hypocenter,  and  the 
resulting  internal  consistency  of  the  fit  to  the  data  may  be  deceptively  good.  For  explosions 
with  known  hypocentral  parameters,  for  example,  it  usually  happens  that  the  true  location  lies 
far  outside  the  estimated  uncertainty  bounds,  indicating  that  such  bounds  may  be  severely 


misleading.  Therefore,  we  have  chosen  to  estimate  the  probable  errors  in  the  data  a  priori,  on 
the  basis  of  experience  with  events  of  known  location,  etc.  These  estimates  will,  of  course,  be 
considerably  larger  than  the  actual  measurement  errors,  or  than  the  standard  deviations  of  the 

data  relative  to  a  best-fit  hypoccnter,  but  they  will  more  realistically  represent  the  true  uncer¬ 
tainty  in  the  location. 

The  basic  method  for  calculating  locations  is  a  simple  extension  of  the  method  first  described 
by  Bolt.  The  (nonlinear)  relation  between  the  hypocentral  parameters  and  each  observed  seismic 
datum  (arrival  time  or  slowness)  is  expressed  in  linearized  form  in  terms  of  changes  from  some 
initial  values.  For  an  arrival  time  datum,  we  get 

6t  =  6to'o-o;60o+^^o  S-^  "«o  (1-1) 

whore 


t  =  arrival  time 


t  =  origin  time 


Oo  =  epi central  colatitude 
</>Q  =  epicentral  longitude 
h  =  focal  depth. 

1  ho  partial  derivatives  in  Fq.  (1-1)  are  given  by 
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where 


T(A,  ho)  =  travel  time 


A  =  epicentral  distance 


tQ  -  azimuth  from  epicenter  to  station. 


A  dt/dA-azinnith  measurement  is  considered  as  a  vector  slowness: 


Ot  .  . 

PX  =  ^A  Slnf 


ot 

Py  =  0  A  COS  t 


where  t  is  the  azimuth  from  the  station  to  the  epicenter.  Corresponding  to  Fq  (1-1),  we  get 
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The  necessary  partial  derivatives  are 
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and  similarly 
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In  the  foregoing,  we  have  made  use  of  the  relations 
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Given  arrival  time  and/or  dt/dA-azimuth  data  for  a  seismic  event  and  a  preliminary  location, 
then  Rqs.  ( I  —  1 )  and  (1-2)  provide  us  with  a  (generally  overdetermined)  system  of  linear  algebraic 
equations  to  be  solved  for  the  desired  changes  in  the  hypocentral  coordinates.  Written  in  matrix 
form 


A  x  =  b 


I 


where-  A  is  a  matrix  of  the  partial  derivatives  appearing  in  Kqs.  (1-1)  and  (1-2)  and  lias  the  form 


1  /a. 
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datum 


A 


0 


0 


slowness 
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Notice  that  we  have  weighted  each  equation  inversely  with  the  estimated  error  a  for  the  corre¬ 
sponding  datum,  x  is  a  column  vectoi  of  1  he  desired  changes  in  the  hypocentral  coordinates: 
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b  is  a  column  vector  of  the  residuals  between  the  observed  data  and  the  values  calculated  for  the 
trial  hypoccntcr  (again  weighted  inversely  by  the  estimated  errors): 

r  6^.1 


h 


6p  fa . 

x. |  /  J 

6  p  fa 

y. i/  j 


7 


The  square'  ol  the  modulus  of  li  is  the  conventional  measure  of  goodness  of  fit  y2; 


\2  =  hTb 


(T  indicates  transposition).  (Ft  should  be  emphasized  that  we  are  taking  care  not  to  interpret 


X2  too  literally  in  terms  of  probabilities.  Standard  tables  of  the  y2  distribution  are  based  upon 


the  assumption  of  independent,  normally  distributed  errors  which,  as  mentioned  above,  does  not 


apply  here.  We  use  y  merely  as  a  rough  estimate  of  whether  the  errors  are  about  the  size  we 


expected  them  to  be.)  Since  the  system  of  Kqs.  (1-3)  is  overdetermined,  it  cannot  lie  solved  ex¬ 
actly  and  we  seek  instead  the  solution  which  minimizes  (Ax  -  b]1’  (Ax  -  b],  which  is  the  solution 


to  the  even-determined  system 


A  F  A  x 


A  rb 


(1-4) 


furthermore,  since  Kqs.  (1-3)  are  merely  linearized  approximations,  they  must  be  solved  itera¬ 
tively,  each  time  with  the  revised  hypocentral  parameters,  until  these  parameters  cease  to  change. 
To  estimate  the  uncertainty  in  the  final  location,  we  need  to  know  the  effect  upon  the  data 


residuals  ol  possible  changes  in  the  hypocentral  parameters.  This  problem  may  be  analyzed 
conveniently  in  terms  of  the  "shifted  eigenvalue  problem"  for  nonsquare  matrices3*4: 


AV  -  A  r 


a1  r  =  A  V 


(1-5) 


where  l  and  V  are  matrices  whose  columns  are  eigenvectors,  and  A  is  a  diagonal  matrix  whose 

elements  are  eigenvalues.  The  eigenvectors  are  guaranteed  to  be  orthogonal  to  each  other,  and 
are  assumed  to  be  normalized  to  unit  length,  from  Kq.  (1-5)  we  see  that 


A  1  AY  =  A2V 


so  that  V  and  A^may  be  calculated  using  conventional  eigenvalue-eigenvector  programs  and  the 
square  matrix  A  A  appearing  in  Kq.  (1-4).  Now,  from  the  orthonormality  of  the  eigenvectors,  it 


follows  that  A  can  lie  factored  as 

yT 


A  =  11 A  V 

if  we  express  a  hypoc enter  perturbation  x  in  terms  of  the  eigenvectors 


x  =  V  y 

then  tlie  corresponding  change  h  in  the  residual  vector  is 


b  =  A  x  =  UA  y 


(1-6) 


This  change  is  orthogonal  to  the  error  vector  b  corresponding  to  the  best-fit  location  since  by 

I<’m  11-41  - „r  .....  _ _ _  ...  -  '  J 


Kq.  d-4),  convergence  of  the  iterative  process  implies  bQ  is  orthogonal  to  the  columns  of  A. 
^  herctoro,  tlic  modulus  ot  1j  ^jivos  tlic  inc rouse  in  ^ 


liquation  (1-6)  says  that  a  change  of  the  hypocentral  parameters  in  the  direction  of  an  eigen¬ 
vector  v.  and  of  length  i/\.  will  increase  y2  by  unity.  Thus,  the  eigenvectors  and  eigenvalues 
completely  specify  the  region  of  the  four-dimensional  space  of  location  and  origin  time  which  will 
satisly  the  data  within  any  desired  accuracy.  Quality  information  of  this  sort  is  indispensable  to 
the  proper  interpretation  of  any  calculated  hypocenter.  Particularly  for  small  events,  the 
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rcsolut  on  of  the  available  data  is  often  so  poor  that  a  calculated  location  without  supplementary 
inform  tion  on  its  uniqueness  may  be  worse  than  useless. 


B.  K,  Julian 


D  FORTRAN  SYSTEM  MODIFICATIONS 


The  Fortran  system  for  the  two  PDP-7  computers  has  been  modified  to  take  advantage  of 
recently  installed  floating-point  hardware.  Most  of  the  programs  used  at  this  facility  show  a 
marked  improvement  in  calculation  speed.  Measurements  taken  on  a  few  programs  with  limited 
input/output  showed  improvement  by  about  a  factor-of-5,  Measurements  on  some  representative 
actual  Fortran  routines  showed:  arithmetic  calculations  5  times  faster,  trigonometric  functions 
about  13  times  faster,  and  exponential  functions  17  times  faster, 

P.  A.  Neilson 
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11.  NORSAR 


A.  SllORT-PERlOD  CHARACTERISTICS  OP  NORSAR 

In  the  last  SATS  (31  December  1912,  DDC  AD-757560),  examples  of  amplitude  variations 
across  NORSAR  from  teleseismic  events  illustrated  the  complexity  of  scattering  produced  by 
deep  structure  under  NORSAR.  The  measured  amplitudes  consisted  of  peak-to-peak  motion  in 
nanos/second  measured  on  the  center  sensor  trace  of  each  .subarray. 

One  of  the  puzzling  features  of  the  data  was  that  inconsistent  patterns  across  NORSAR  were 
obtained  for  two  closely  located  presumed  explosions  in  Eastern  Kazakh.  This  result  implied 
that  the  amplitude  patterns  were  controlled  not  by  scattering  of  waves  by  structure  under  NORSAR, 
but  rather  by  source  and  transmission  path  effects  further  removed  from  NORSAR.  This  was 
contrary  to  results  obtained  for  two  deep  Hindu  Kush  events  with  close  epicenters.  Although  the 
NOS  dep.hs  for  the  Hindu  Kush  events  difter  by  80  km,  very  similar  amplitude  patterns  were 
produced  across  NORSAR.  These  data  suggested  that  the  deep  structure  under  NORSAR  was  the 
dominant  factor  which  controls  the  amplitude  pattern  across  the  array. 

These  contradictory  results  have  been  resolved  by  a  more-detailed  frecpiency  domain  analysis 
of  the  data.  The  main  conclusion  is  that  the  amplitude  patterns  across  NORSAR  are  strongly 
1  requeney  dependent,  and  that  measurements  of  nanos/second  from  records  give  consistent  ampli¬ 
tude  patterns  only  for  source  functions  with  similar  spectra  and  phase.  Events  from  the  same 
regior.r  with  dissimilar  spectral  peaks  may  have  very  different  amplitude  patterns  in  nanos/second 
across  NORSAR.  This  scattering  effect  is  accentuated  at  NORSAR  by  the  crustal  geology.  The 
shield  under  NORSAR  has  a  very  high  Q  which  transmits  a  wide  band  of  frequencies  to  the  sensors. 


TABLE  ll-l 

NOS  DATA  FOR  TWO  PRESUMED  SOVIET  EXPLOSIONS 


Event 

Date 

(1971) 

Origin  Time 
(GMT) 

Location 

Magnitude 

(mb) 

1 

6  June 

04:02:57.  1 

50. 0°N  77. 8°E 

5.5 

2 

30  June 

03:56:57.2 

50. 0°N  79.  1°E 

5.4 

As  an  example,  we  present  here  the  amplitude  patterns  of  two  presumed  explosions  in 
Eastein  Kazakh  which  have  clearly  unequal  source  functions.  The  NOS  data  for  the  events  are 
given  in  Table  11-1 . 

Figures  11-1  and  11-2  show  the  center  subarray  sensors  and  steered  beams  at  NORSAR  for 
events  1  and  2,  respectively.  Event  1  was  chosen  because  it  has  an  impulsive  high-frequency 
beam  with  a  period  of  ~0,4  sec.  Event  2  has  a  more-complex  beam  with  an  apparent  period  of 
~  1  -0  sec.  I  he  spectra  of  these  two  beams  are  displayed  in  Fig.  11-3,  The  beam  spectrum  of 
event  1  is  nearly  flat  from  1  to  3  Hz,  whereas  the  spectrum  of  event  2  peaks  at  1  llz  with  a  pro¬ 
nounced  notch  at  ~1.7  Hz,  due  to  the  interference  seen  on  the  beam  in  Fig.  11-2, 

Evert  1  can  be  used  to  calibrate  the  array  if  we  assume  that  the  beam  is  a  good  estimate  of 
the  source  incident  to  the  crustal  structure  under  NORSAR.  The  spectrum  of  each  sensor  wave¬ 
form  for  this  event  will  indicate  the  spectral  response  of  the  structure  under  NORSAR  to  Eastern 
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Kazakh  events.  In  Fig.  11-4,  the  spectral  responses  of  the  individual  waveforms  of  event  1  are 
plotted  at  frequencies  of  1,  2,  and  3  llz.  These  speetra  were  eomp  .ited  using  20  sec  of  eosine- 
tapered  data.  On  the  right-hand  side  of  the  figure  are  the  spectral  values  of  the  steered  beam. 

Figure  11-4  clearly  shows  that  the  array  response  varies  from  1  to  2  llz.  Sensors  16,  17 
and  18,  for  example,  have  nearly  equal  responses  at  1  Hz,  whereas  at  2  Hz,  sensor  17  has  a 
response  well  below  that  of  16  or  1  8.  This  explains  the  apparent  inconsistency  in  waveform 
shapes  at  these  sensors  for  events  1  and  2.  In  Fig.  11-1,  the  waveform  amplitudes  at  sensors 
1 6,  17,  and  18  match  the  2 -llz  pattern  at  the  same  sensors  in  Fig.  11-4.  However,  in  Fig.  11-2 
the  same  sensors  have  nearly  equal  amplitudes,  and  match  the  1-Hz  pattern  in  Fig.  11-4.  Even 
the  high  frequency  content  of  sensors  16  and  18,  which  is  absent  on  sensor  17  in  Fig.  11-2,  is 
predicted  by  the  array  response  in  Fig.  11-4. 

In  order  to  investigate  the  consistency  of  crustal  scattering  in  more  detail,  least-squares 
time-domain  filters  were  computed  for  shaping  the  beam  of  each  event  to  the  individual  sensor 
waveform.  This  removes  the  source  from  each  waveform  and  isolates  the  scattering  received 
at  each  sensor.  These  filters  contain  smoothed  phase  and  amplitude  information  in  the  frequency 
domain.  Figures  11-5  and  11-6  show  the  computed  filters  for  events  1  and  2,  respectively.  In 
spite  of  the  different  source  spectra  of  the  two  events,  the  two  sets  of  filters  are  very  similar 
in  shape  and  amplitude. 

In  Fig.  11-7,  some  of  the  more  coherent  pairs  of  filters  for  the  two  events  are  shown  aligned 
in  time.  The  filters  appear  to  be  most  coherent  in  the  first  one  or  two  seconds  of  time,  although 
in  some  cases  coherent  bursts  arrive  later  in  time.  Coherence  in  the  time  domain  carries  over 
into  the  frequency  domain.  In  Fig.  11-8,  the  speetra  of  the  two  filters  for  sensor  2  are  superposed. 
The  notch  in  amplitude  at  2  Hz  in  both  speetra  indicates  that  the  scattering  recorded  at  this  sensor 
is  probably  a  multipath  phenomenon. 

Attempts  are  being  made  to  examine  the  spectral  amplitudes  at  NORSAR  from  events  of 
various  azimuth  and  distance  and  relate  these  to  erustal  and  Moho  structure. 

C.  W.  Frasier 

B.  NORSAR  ARRAY  MISLOCATIONS 

A  sizable  data  base  of  good-quality  well-recorded  events  at  NORSAR  is  being  used  to  obtain 
values  for  the  array  epicenter  location  errors,  as  well  as  the  subarray  time  anomalies.  While 
the  complete  pieture  of  the  time  anomalies  is  not  yet  available,  we  have  obtained  a  good  pieture 
of  the  epicenter  misloeations  that  NORSAR  would  produce  without  using  any  station  corrections. 
The  time  delays  were  obtained  from  the  event  processor  at  NORSAR  and  were  corrected  to  agree 
with  visual  recordings.  A  plane  wave  was  fitted,  in  a  least-squares  sense,  to  eaeh  set  of  delay 
times.  The  resulting  measurements  of  velocity  and  azimuth,  and  therefore  latitude  and  longitude, 
were  plotted  against  the  true  epicenter  or  velocity  and  azimuth  obtained  from  the  standard  J-B 
tables.  The  resulting  plots  show  effects  that  were  not  widely  observed  at  LASA,  yet  seem  to  be 
prevalent  at  NORSAR. 

On  first  analysis  of  the  array  diagram  in  Fig.  11-9  (a  polar  plot  with  dT/dA  as  the  radial 
component),  there  is  no  one  overriding  trend  as  was  observed  at  LASA.  In  this  figure,  the  slow¬ 
ness  and  azimuth  of  the  true  epicenters  reported  by  NOS  are  represented  by  the  arrowheads  and 
plotted  at  the  slowness  values  given  from  the  J-B  tables.  The  array  measurements  of  azimuth 
and  slowness  are  plotted  as  the  tails  of  the  arrows.  Only  P-wave  data  were  used.  It  is  difficult 
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to  see  any  universal  bias  to  the  arrows  in  Fig.  II-9.  This  was  not  the  ease  at  LASA,1  where  there 
appeared  to  be  a  constant  misloeation  term.  Figure  11-10  is  a  projection  of  the  coastal  outlines 
into  slowness  and  azimuth  and  will  help  in  interpreting  the  array  diagram.  On  a  large  projection 
of  the  array  diagram,  it  is  evident  that  there  are  a  number  of  regions  where  the  misloeation 
vector  will  ehange  rapidly  over  small  ehanges  in  azimuth  or  slowness.  These  areas  are  not 
usually  the  same  areas  of  rapid  ehange  observed  at  LASA. 

When  the  misloeation  vectors  are  plotted  on  a  geographic  map,  it  is  possible  to  see  the  areas 
of  rapid  ehange  more  clearly.  Figure  11-11.  for  example,  shows  how  the  loeatic  n  error  ehanges 
as  the  epicenters  move  from  western  United  States  to  eentral  South  America.  Oier  this  region, 
the  azimuthal  constituent  of  the  location  error  varies  rapidly  while  the  slowness  values  of  these 
events  decrease  right  into  the  core.  What  we  are  seeing  here  is  an  array  location  error  changing 
sign  and  passing  through  a  zero-loeation  error  point  around  Central  America.  This  is  quite  clear 
when  viewed  in  Fig.  II-9  at  the  azimuth  of  270".  A  small  set  of  events  in  the  Central  America 
region  shows  still  another  interesting  effect.  These  events  from  the  Oaxaca  district  of  Mexieo, 
indicated  by  the  box  in  Fig.  H-ll,  are  actually  quite  close  to  each  other  and  yet  with  only  a  slight 
shift  in  epicenter,  less  than  7  5  km,  the  array  locations  are  shifted  by  as  much  as  700  km.  One 
can  draw  a  llne  through  this  region  which  will  clearly  separate  the  epicenters  into  two  distinct 
anomalies,  yet  each  epicenter  will  be  within  a  few  tens  of  kilometers  of  the  line.  This  feature 
clearly  seems  to  be  source  or  path  related,  whereas  the  overall  trend  within  Fig.  11-11  could 
possibly  be  due  to  near  array  effects. 

Another  interesting  area  to  note  is  in  western  Asia  around  the  China,  Russia,  Afghanistan 
borders  (Fig.  11-12).  In  this  region,  where  the  seismic  zones  of  Hindu  Kush  and  northern  Sinkiang 
meet,  the  azimuth  from  NORSAR  will  change  by  only  a  few  degrees,  yet  the  epicenter  location 
error  will  almost  completely  reverse  itself.  It  is  supposed  that  this  feature  is  also  primarily  a 
souree/path  effect. 

As  more  processing  of  the  data  occurs,  it  is  expected  that  many  more  interesting  features 
will  emerge.  Still  to  be  analyzed  in  detail  however  are  the  individual  subarray  anomalies. 


R. M.  Sheppard 


C.  ANALYSIS  OF  40-SEC-PERIOD  RAYLEIGH- WAVE  CODA  AT  NORSAR 


The  occasional  long  persistence  of  Rayleigh-wave  coda  can  make  it  difficult  to  apply  the 
Ms’mb  discriminant  to  subsequent  signals  received  at  a  seismic  station.  Previous  work2’3  at 
LASA  revealed  the  contribution  of  multipath  energy  to  this  coda  and  the  desirability  of  studying 
*40  “see  energy  which  should  suffer  less  multipathing. 

A  similar  analysis  of  the  40-see-period  Rayleigh-wave  coda  is  being  done  at  NORSAR.  Some 

_ • _ ....  ..  -i  . 


preliminary  results  are  now  presented.  A  total  of  16  events,  each  with  m,  *  5.5,  located  at 


various  azimuths  and  distances  from  NORSAR,  have  been  selected.  Therewere  no  events  occur¬ 
ring  close  in  time  to  the  chosen  events,  according  to  the  NOS  monthly  summary  reports  and  the 
NORSAR  bulletin.  The  high-resolution  wavenumber  analysis  method  was  modified  so  as  to  pre¬ 
filter  the  NORSAR  L.PZ  seismometer  data  in  a  manner  similar  to  that  done  previously  at  LASA.2 

I  he  rcsilIts  of  thc  analysis  for  the  12  .lime  1972  Aleutian  Island  event  are  shown  in  Figs  U-l  3(a-d) 
and  II-14(a-d). 

In  Fig.  11-1  3(a)  we  see  that,  during  the  first  200  see  following  the  onset  time  of  the  40-sec 
period  group,  the  peak  power  occurs  at  an  azimuth  of  0"  and  with  a  phase  velocity  corresponding 
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approximately  to  that  of  a  fundamental-mode  Rayleigh  wave.  Thus,  during  this  time,  the  group 
arrives  at  NORSAR  from  very  close  to  the  true  azimuth  of  this  event  at  NORSAR,  which  is  358.5°. 
In  Figs.  II— 1  3(b)  and  (c),  there  are  very  clear  indications  of  multipath  propagation  of  40-sec  period 
groups  at  NORSAR.  An  effort  is  under  way  to  try  to  identify  the  propagation  paths  for  the  various 
phases  which  are  being  detected  in  the  coda  at  NORSAR, 

Preliminary  results  at  NORSAR  seem  to  indicate  that  the  coda  power  level  decreases  more 
rapidly  than  that  found  at  LASA.  For  example,  the  data  in  Figs.  II-l  3  and  11-14  show  that  the 
power  level  decreases  by  30  dB  after  20  minutes.  This  is  about  10  dB  better  than  the  average 
value  found  at  LASA,2  which  indicates  that  NORSAR  may  have  an  important  role  to  play  in  a  net¬ 
work  for  detecting  interfering  Rayleigh  waves. 

J.  Capon 
H.  Bungum* 
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Fig.  II- 9.  NORSAR  array  diagram  for  P-waves. 
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Fig.  11-10.  Coastal  outlines  in  slowness  projection  from  NORSAR. 


Fig.  11-11.  NORSAR  location 
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111.  LATERAL  IlETEllOC.ENEIT V 


A.  PREDICTING  RAVLElGll- WAV E-PROPAC.ATION  BEHAVIOR 
FROM  GEOLOGICAL  AND  GEOPHYSICAL  DATA 

In  the  SATS  dated  30  June  1972  (DDC  AD-748304),  preliminary  resuLs  were  given  of  an 
experiment  aimed  at  using  ray  tracing  to  explain  in  detail  multipath  propagation  and  lateral  re¬ 
traction  of  Rayleigh  waves.  The  results  of  that  experiment,  based  on  a  very  crude  model  con¬ 
taining  only  oceans  and  continents,  compared  well  with  observed  data  and  have  encouraged  us  to 
undertake  experiments  with  more  complicated  models  to  find  out  the  degree  of  detail  with  which 
surface  wave  propagation  on  the  real  earth  may  be  predicted. 


There  are  several  ways  in  which  a  detailed  understanding  of  surface  wave  propagation  would 
be  ol  value  in  the  seismic  discrimination  problem.  Surface  wave  magnitude  (M  )  measurements 
are  allected  bv  regional  differences  in  dispersion  and,  probably  more  severely,  by  geometric 
spreading  caused  by  lateral  refraction.  Accurate  path  corrections  for  these  effects  would  make 
Ms  measi,roments  lor  small  events  more  reliable.  The  same  dispersion  and  spreading  informa¬ 
tion  would  enable  synthetic  seismograms  to  be  computed  much  more  accurately  than  is  now  pos¬ 
sible,  and  would  thus  improve  our  ability  to  detect  small  signals  hy  matched  filtering.  Multi- 
pathing  calculations  would  enable  reliable  associations  to  be  made  between  observed  signals  and 
seismic  events  in  cases  of  ambiguity,  and  thus  are  of  great  potential  use  in  the  interfering  event 
problem. 


We  have  chosen  to  deal  with  the  western  United  States,  because  relevant  surface  wave  data 
are  available  from  l.ASA  and  because  extensive  geological  and  geophysical  data  on  crust  and 
upper  mantle  structure  are  available.  The  results  discussed  here  apply  to  models  incorporating 
inlormation  on  crustal  thickness  (see  Fig.  1II-1),  upper  mantle  velocities,  and  crustal  structure 
differences  between  geological  provinces. 

Figures  lll-2(a)  through  (d)  show  Rayleigh  wave  rays  traced  outward  from  LASA  in  such  a 
model.  These  calculations  may  be  compared  with  the  observations  reported  by  Landers  in  the 
last  SATS  (Fig.  1-19)  of  azimuth  anomalies  at  LASA  for  Rayleigh  waves  from  the  Nevada  Test 
Site.  At  10 -sec  periods,  the  waves  arrived  10°  to  the  west  of  the  great  circle  azimuth,  in  nearly 
perfect  agreement  with  the  prediction  of  Fig.  Ill-2(a).  This  anomaly  is  caused  by  refraction  at 
the  boundary  ol  the  Snake  River  Plain,  in  which  the  crust  (consisting  entirely  of  basalt)  has  very 
high  seismic  velocities.  The  observed  azimuth  anomaly  disappears  at  about  20  sec,  also  in 
agreement  with  prediction  [Fig.  Ill— 2(b)  ],  but  at  longer  periods  the  observed  waves  arrive  to  the 
east  of  the  great  circle,  and  no  such  effect  is  predicted  [Figs.  III-2(c)  and  (d)  J.  An  attempt  to 
resolve  this  discrepancy  by  raising  the  mantle  velocity  in  the  Colorado  Plateau  proved 
unsuccessful. 

Thus,  on  the  basis  of  these  limited  data,  we  have  predicted  the  short-period  behavior  of 
Rayleigh  waves  quite  well,  and  discovered  that  our  knowledge  of  the  Colorado  Plateau-Basin  and 
Range  boundary  is  inadequate  in  some  yet  unknown  respect.  Current  work  is  aimed  at  refining 
the  model  (e.g.,  incorporating  sediment  thicknesses)  and  acquiring  more  data  against  which  to 
test  it. 

B.  R.  Julian 

B.  BODY-WAVE  MULTIPATH1NG  -  EVIDENCE  FROM  TIME-DOMAIN  ANALYSIS 
\ 

In  the  last  SATS,  Capon  and  I  presented  evidence  which  indicated  that  waves  from  a  region 
in  South  America  followed  paths  from  source  to  receiver  which  were  dependent  on  frequency. 
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We  four  J  in  particular  that  frequencies  less  than  1  Hz  followed  a  path  significantly  different  in 
azimuth  from  that  followed  by  frequencies  greater  than  1  Hz. 

A  very  simple  and  large  earthquake  in  the  Hindu  Kush  (20  January  1972,  depth  =  220  km, 
mfa  =  6.0)  has  revealed  multipathing  whieh  is  visible  even  on  the  unprocessed  seismometer  out¬ 
puts.  Figure  111-3  shows  6  subarray  sums  for  short-period  P-waves  from  aeross  the  whole 
aperture  of  the  array,  It  is  clear  that  the  signal  consists  of  a  primary  pulse,  followed  by  later 
aetivity  among  which  a  secondary  pulse  is  prominent  and  has  been  aligned  in  the  figure.  Align¬ 
ment  of  the  secondary  pulse  however  yields,  as  is  obvious  from  Fig.  111-3,  a  different  arrival 
vector  from  alignment  of  the  primary  pulse. 


The  two  phases  elearly  arrive  with  different  values  of  dT/dA,  although  the  azimuths  along 
which  they  arrive  arc  not  significantly  different.  Note  that  the  determinations  of  wave  parameters 
for  Table  111-1  were  made  by  least-squares  fitting  of  a  plane  wave  to  the  raw  arrival  times.  No 
subarray  corrections  were  applied.  For  this  reason,  although  the  data  in  Table  111-1  represent 
estimates  relatively  free  of  personal  prejudice,  the  absolute  values  of  dT/dA  and  azimuth  should 
not  be  taken  seriously,  as  they  are  sub  ject  to  bias  from  structure  beneath  the  array. 

It  is  clear  that  the  dT/dA  values  are  too  significantly  different  for  us  to  consider  the  possi¬ 
bility  that  this  is  a  source  phenomenon,  such  as  a  double  rupture.  Examination  of  the  single- 
pulsed  PKKP  at  L.ASA  confirms  this.  Nor  can  it  be  easily  ascribed  to  the  region  under  LASA, 

If,  for  instance,  it  were  a  P-to-S  conversion  phenomenon,1  we  would  expect  that  the  delays  be¬ 
tween  phases  1  and  2,  representing  sediment  thicknesses,  would  be  random,  This  they  are  not, 

In  Fig,  111-4  we  plot  time  residuals  to  phase  2  against  time  residuals  to  phase  1  for  each  subarray, 
and  it  is  clear  that  they  are  strongly  correlated.  This  implies  that  the  two  waves  have  experienced 
the  same  perturbations  on  a  plane  wave  front  on  approaching  LASA,  and  therefore  that  the  effect 
is  not  locally  generated. 

Further  powerful  constraints  on  the  location  of  the  source  of  multipathing  eome  from  the 
time  delay  of  2  sec  between  the  first  and  second  pulse.  Figure  111-5  is  a  very  schematic  ray 
diagram.  Returning  from  LASA  are  a  pair  of  rays  separated  by  an  (exaggerated)  angle  intended 
to  represent  the  dT/dA  difference  between  the  two  pulses.  Note  that  the  second  pulse  has  a 
larger  dT/dA,  bottoms  at  a  shallower  depth,  and  therefore  cannot  be  PcP.  If  the  earth  were 
laterally  homogeneous  and  a  sharp  turn  were  introduced  into  ray  paths  only  by  a  localized  feature 
which  did  not  affect  the  travel  time  along  the  path,  to  produce  a  2-sec  time  differential  it  would 
be  necessary  for  the  two  ray  paths  to  run  back  until,  within  600  km  of  the  surface  near  source, 
one  of  the  rays  turned,  This  is  one  of  the  options  sketched  in  Fig,  111-5;  another  is  that  there 
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are  sufficiently  strong  velocity  gradients  just  above  the  core  that  there  is  triplication  in  the  travel 
time  curve.  A  variety  of  other  explanations  intermediate  between  these  two  extremes  exists,  but 
it  should  be  clear  that  any  complex  structure  beneath  LASA  (say,  on  the  650-km  discontinuity), 
although  it  might  be  capable  of  generating  multipaths,  would  not  in  any  reasonable  model  yield 
time  delays  of  2  sec  or  more.  I  thus  conclude  that  a  complexity  of  earth  structure  exists  back 
along  tbe  ray  paths  either  in  the  deep  mantle  or  beyond  toward  the  source,  such  that  at  least  two 
possible  paths  exist  from  source  to  receiver. 

D.  Davies 


C.  AMPLITUDES  OF  P-WAVES  THAT  HAVE  TRAVERSED  THE  DEEP  MANTLE 

We  have  studied  the  amplitudes  of  short-period  P-waves  in  the  distance  range  90°  to  105° 

and  some  rather  striking  features  have  emerged.  A  full  account  of  this  \  ork  is  given  in  a  paper 
2 

now  in  press,  so  we  will  simply  outline  the  main  conclusions. 

Seismicity  around  the  Pacific  from  New  Guinea  to  Kermadec  gives  a  rather  good  coverage 
at  the  distance  98°  to  102°  for  azimuths  from  235°  to  275°  from  LASA.  it  is  clear  that  events 
of  the  same  body  wave  magnitude  (NOS)  do  not  necessarily  have  the  same  amplitude  on  an  equi¬ 
distant  contour.  For  instance,  "visibility"  to  Santa  Cruz  is  much  inferior  to  visibility  to  points 
in  Fiji  and  the  Solomon  Islands  at  the  same  distance.  It  is  possible  to  eliminate  effects  under 
LASA  fairly  easily,  but  more  attention  must  be  paid  to  the  possibility  of  source  region  effects. 
However,  after  careful  consideration  of  the  effect  of  a  dipping  slab,  or  a  low  Q  zone,  we  have 
concluded  that  the  source  region  cannot  be  the  culprit  either.  We  are  thus  left  with  propagation 
effects  in  the  middle  or  deep  mantle  as  the  most  likely  cause.  It  is  not  immediately  possible  to 
distinguish  between  velocity  and  Q  heterogeneity.  The  dimension  of  the  anomaly  at  LASA  is  10° 
of  azimuth.  This  suggests  some  feature  a  few  hundred  kilometers  across  deep  in  the  earth. 

R.  E.  Needham 
D.  Davies 

D.  ANALYSIS  OF  ARRAY  DIAGRAM  AT  LASA 

■[ 

In  the  last  SATS,  a  discussion  was  presented  on  the  representation  of  approach  angle,  or 
slowness,  anomalies  for  P-wave  signals  at  an  array  in  slowness  (dT/dA,  azimuth)  space.  Such 
a  representation,  known  as  an  array  diagram,  involves  the  plotting  of  these  anomalies  as  arrows. 
The  tail  represents  the  dT/dA  and  azimuth  measured  for  the  P-wave  signal  at  the  array,  and  the 
head  represents  the  dT/dA  and  azimuth  that  would  be  obtained  from  the  true  (or,  in  this  case, 

NOS)  location  and  ,J-B  tables.  This  approach  has  been  used  by  Manchee  and  Weichert,3  Vinnik 

4  5 

and  Nikolayev,  and  by  Davies  and  Sheppard  at  LASA. 

If  the  anomaly  were  caused  by  a  plane-dipping  interface  directly  under  the  array,  the  array 

diagram  would  consist  of  a  set  of  error  vectors  whose  magnitudes  were  all  approximately  the 

same  and  which  all  pointed  in  the  direction  of  dip  of  the  interface.  Thus,  the  effect  of  a  sloping 

interface  can  be  obtained  from  the  array  diagram  by  fitting  a  constant  vector  to  the  error  vectors 

3 

in  the  diagram.  These  results  were  first  pointed  out  by  Manchee  and  Weichert.  Davies  and 
Sheppard5  pointed  out  that  a  smoothly  varying  distribution  of  arrows,  all  approximately  of  equal 
length  and  pointing  in  the  same  direction,  would  indicate  a  bias  introduced  beneath  the  array  - 
not  necessarily  from  a  dipping  interface  out,  for  example,  from  a  distribution  or  coating  of 
delays  that  are  not  directionally  independent. 

The  spatial  rate  of  variation  of  the  error  vectors  can  at  least  formally  be  associated  with 
heterogeneities  at  certain  depths  in  the  earth.  The  more  rapid  these  spatial  variations  become, 
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e  more  distant  from  LASA  the  heterogeneities  must  be  placed  to  explain  the  anomalies  The 
reason  for  this  can  be  seen,  qualitatively  at  least,  on  the  basis  of  the  separation  of  the  ray  paths 
rnm  the  various  event  epicenters  to  LASA.  We  now  wish  to  give  a  more  quantitative  proof  for 
is  important  result,  since  it  forms  the  basis  for  an  analysis  of  the  array  diagram.  Suppose 
•  at  the  region  under  LASA  can  be  divided  into  a  series  of  sub-regions,  each  of  which  can  be 
considered  a  horizontal  disk  extending  between  certain  depths  as  shown  in  Fig.  HI-6,  The  outer 
boundaries  ot  these  disks  are  defined  by  the  P-wave  ray  paths  which  have  a  dT/dA  of  abou 
10  sec/deg,  as  indicated  in  Fig.  HI-6.  We  consider  first  the  topmost  disk  extending  betwe  0 
and  75  km  The  effect  of  this  region  on  the  propagation  of  the  P-waves  can  be  obtained  by  con¬ 
sidering  the  region  to  be  approximated  by  a  combination  of  a  time  delay  coating,  an.  a  randjri 
medium,  it  was  shown  in  the  last  SATS*  that  this  approximation  appears  to  be  reasonably  good 
or  the  earth  structure  under  LASA.  In  addition,  it  was  shown  that  LASA  tends  to  average  out 
the  effects  due  to  the  laterally  heterogeneous  medium  on  the  P-wavv  *ig  ,1„,  ^  tfmt  lL 

foVap0031!"0  dC,nreS-  ThlS  COating  Produces  an  anomaly  component  which  must  be  the  same 
all  vectors  in  the  array  diagram,  since  it  is  located  so  close  beneath  the  array  that  all  P-wavt 

I  ay  paths  must  pass  through  it.  Therefore,  the  effect  of  the  coating  may  be  obtained  by  deter-  '  ' 

Z  Fie  ill  °f  aU  thG  err°r  VeCt°rS  ^  the  array  diaRram>  K  mUSt  b°  said  as  a  caution 

that  Fig.  IH-6,  showing  rays  as  straight  lines,  is  clearly  no  more  than  a  schematic  representation 
intended  to  start  discussion.  P  n  t 

We  now  consider  the,  ,h,  erred  .he  region  i„  the  0-  ,5-ltm  „nge  has  been  removed 

may  diagram  by  subtracting  the  average  ot  all  the  error  vectors  trom  each  or  those  vectors 
We  now  examine  tin,  j(»k  p  ih*  7^-  ^  or  iS[  ■  „  h 

„  «t.  ra„fee.  1  ne  nonzontal  cross  section  of  this 

chsk,  at  a  depth  of  about  1  50  km.  is  equal  to  about  four  times  the  area  of  the  LASA  aperture 

18  n0WaSSUm°d  that  thiS  **  be  divid*d  foul  regions,  each  of  wmch  can  once  again  be 
approximated  by  a  tandem  combination  of  a  coating  and  random  medium.  Since  there  are  four 
such  coatings,  there  will  he  four  independent  error  vectors  introduced  in  the  array  diagram. 

e  .  ,  .  ,  -t,,  ,  hOtors  iSCarTSed  by  an  associated  coating,  since  the  effect  of  each  of  the 
random  media  is  averaged  out  by  LASA  as  indicated  previously.  In  addition,  each  vector  would 
e  approximately  constant  over  an  area  roughly  equal  to  one-fourth  of  the  total  area  of  the  array 
diagram^  Thus,  in  order  to  evaluate  each  of  these  vectors,  the  array  diagram  is  divided  into 
our  quadrants  and  an  average  of  the  error  vectors  (with  the  average  vector  removed)  is  per- 
orme  ln  eac  quadrant.  Once  again,  we  can  consider  that  the  effect  of  the  region  in  the  75-  to 

from Th  ranSe  T  been  rem°VCd  fr°m  the  3rray  dl3gram  by  subtractin6  each  of  these  four  vectors 

the  arT0dma  7/  ^  ^  ^  rem°Ved’  in  the  ^prUte  quadrants 

of  the  array  diagram.  After  this  has  been  done,  the  effect  of  the  region  in  the  225-  to  375-km 

ePth  range  can  be  considered  and  analyzed  in  a  similar  manner.  The  analysis  can  then  be  ex- 

d /o  ^e  other  regions  shown  in  Fig.  UI-6  extending  to  a  depth  of  975  km.  Thus,  it  might  be 

possible  to  use  this  analysis  to  determine  the  effect  on  the  array  diagram  of  the  crust,  upper- 

and  mid-mantie  structure  under  LASA.  It  should  be  remarked  that  this  is  a  first  attempt  at 

analysis  of  the  array  diagram,  and  is  certainly  not  the  best  possible.  For  instance,  the  array 

diagram  is  expressed  in  terms  of  nonorthogonal  functions  (namely,  means  valid  for  the  whole 

segment).  Nevertheless,  these  early  results  have  some  interest.' 

The  tmelysie  procedure  deserved  previously  has  been  applied  to  the  array  diagram  a.  LASA 

original  data  constated  ot  2253  error  vectors.  However,  in  order  to  provide  uniform  coverage 


in  the  array  diagram,  a  data  decimation  scheme  was  used  in  which  the  array  diagram  was  divided 
into  a  grid  of  20  x  20  cells,  where  each  cell  is  a  square  of  1  x  1  sec/deg.  The  first  ten  vectors 
to  fall  within  a  given  cell  were  taken  as  acceptable  data.  This  led  to  a  total  of  299  error  vectors 
as  shown  in  Fig.  111-7.  The  average  of  these  error  vectors  was  computed  and  found  to  point  in 
the  direction  corresponding  to  an  azimuth  of  15°,  with  a  modulus  or  0.32  sec/deg,  which  is  in 
reasonable  agreement  with  the  corresponding  results  of  Davies  and  Sheppard.5  This  mean  vector 
could  be  associated  with  a  plane-dipping  interface  which  dips  down  at  an  azimuth  of  15°.  These 
parameters  could  be  associated  with  a  mean  dip  of  4°  of  the  Moho  under  LASA,  although  the 
Moho  is  not  necessarily  the  only  source  of  anomaly. 

In  Fig.  111-8,  we  have  the  sum  of  the  average  vector  and  all  the  vectors  obtained  using  the 
sub-averaging  method  described  previously.  The  great  similarity  of  the  array  diagram  in 
big.  111-8  with  that  in  big.  HI -7  should  be  noted.  This  similarity  is  also  borne  out  by  the  small¬ 
ness  ol  the  vectors  sho  'n  in  Fig.IlI-9,  which  depicts  the  difference  between  the  error  vectors 
in  1' igs.  111-7  and  111-8.  The  mean-square  value  of  the  error  vectors  shown  in  Fig.  111-9,  relative 
to  that  in  the  original  array  diagram  shown  in  Fig.  111-7,  is  about  -10.7  dU.  This  means  that  the 
analysis  procedure  has  been  able  to  remove  more  than  90  percent  of  the  original  error  vectors, 
when  measured  in  terms  of  mean-square  value. 

It  is  also  possible  to  find  I  he  mean-square  value  of  the  error  vectors  in  each  depth  range, 
relative  to  that  in  the  original  array  diagram  shown  in  Fig.  111-7.  The  results  of  such  a  computa¬ 
tion  are  shown  in  Fig.  111-10,  where  the  ordinate  has  been  labeled  as  relative  amount  of  hetero¬ 
geneity.  It  is  seen  trom  this  figure  that  there  is  an  approximately  decreasing  amount  of  hetero¬ 
geneity  as  the  depth  increases.  I  his  agrees  roughly  with  our  notion  that  the  earth  becomes  more 
homogeneous  as  the  depth  increases.  A  similar  analysis  has  been  performed  at  NORSAR  and 


the  results  will  be  presented. 


.1.  Capon 
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Fig.  Ill  - 1 .  Elevation  of  Mohorovicic  discontinuity  (km  below  sea  level) 
in  western  United  States,  based  on  all  available  data  (primarily  from 
seismic  refraction  measurements).  Contour  interval  is  5  km. 
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Fig.  HI-5.  Exaggerated  sketch  of  possible  ray  paths  to  generate  multiple 
signals  at  LASA. 
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Fig.  III-6.  Decomposition  of  region  directly  under  LASA. 
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Fig.  Ill— 7.  Original  data  used  in  analysis  of  array  diagram. 
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Fig.  III-8.  Array  diagram  obtained  by  adding  average  and  sub-averaee 
error  vectors. 


IV.  EVASION  AND  COUNTER  EVASION 


A.  MULTIPLE  EXPLOSIONS  THAT  PRODUCE 
EARTHQUAKE- LIKE  SEISMOGRAMS 

Two  short-period  characteristics  of  a  large  class  of  earthquakes  immediately  allow  them  to 
be  discriminated  from  single  explosions.  When  the  rupture  time  of  an  earthquake  exceeds  more 
than  several  seconds  and  when  seismic  energy  is  continuously  released  over  that  time  span, 
seismograms  written  at  teleseismic  distances  will  be  far  more  complex  than  those  from  single 
explosions  even  when  a  large  degree  of  scattering  along  the  explosion-receiver  path  is  present. 
Further,  when  the  depth  phases  pP  and  sP  can  be  identified  and  when  their  delaytimes  indicate 
a  hypocentral  depth  greater  than  is  technically  possible  to  conduct  an  explosion,  then  these  events 
are  also  necessarily  earthquakes  rather  than  single  explosions.  The  possibility  that  these  earth¬ 
quake  characteristics  could  be  successfully  simulated  by  spatially  and  temporally  distributing 
several  explosions  is  investigated. 

A  constraint  on  simulating  complexity  is  that  the  surface-  to  body-wave  energy  ratio  he  con¬ 
sistent  with  the  ratio  that  would  be  obtained  for  an  earthquake.  Thus,  it  is  necessary  to  time 
the  explosions  so  that  the  short-period  body  waves  do  not  combine  constructively,  while  the  long- 
period  surface  waves  do.  Figure  IV- 1  shows  eight  identical  source  functions  with  amplitude 
ratios  of  2,  5,  9,  20,  20,  40.  30,  25  and  time  delays  H  j  1  2,  .3j  <,  7  ,  seC;  alan  aftflWTl 

is  the  sum  of  these  eight  functions.  The  source  function  is  that  of  a  presumed  explosion  from 
Eastern  Kazakh  of  body-wave  magnitude  5.5  as  recorded  at  LASA  (A  =  83”).  The  sum  trace  is 
sufficiently  complex  that  such  a  seismogram  could  not  have  been  caused  by  a  single  explosion. 

The  body-wave  magnitude,  as  conventionally  measured  in  the  4 vyvW,  sv&Hdr 
than  that  for  the  largest  single  explosion  by  0.6.  At  small  epicentral  distances,  surface  waves 
contain  more  high-frequency  energy  and  interference  should  show  up  more  clearly.  Consequently, 
to  examine  the  effect  on  surface  waves,  the  source  Rayleigh  wave  of  an  NTS  explosion  recorded 
at  LASA  is  used  (A  =  12").  Figure  IV-2  shows  the  result  for  the  Rayleigh  wave  with  the  same  de¬ 
lays  and  amplitude  ratios.  The  individual  waves  sum  without  significant  interference,  trebling 
the  maximum  amplitude.  The  result  is  an  increase  in  surface-wave  magnitude  of  0.4.  Fig¬ 
ure  IV- 3  shows  the  total  effect  on  the  Mg:mb  diagram.  The  trend  lines  are  those  determined  by 
Lacoss.  If  the  Mg:mb  value  of  the  largest  single  explosion  occurred  at  A,  then  the  decrease  in 

mL  rnoves  the  point  to  R  and  the  increase  in  Mg  moves  it  to  C,  wfctali  r liLUrJy  Um*  wilhin  tire 
earthquake  population. 

It  is  reasonable  to  conjecture  that  the  longer  the  short-period  waveform  of  a  given  event  the 
greater  the  source  dimensions  of  thot  oy.-nt  aid  cc-nfe quentty,  the  greater  the  Tati,,  ol  Lmg-peTiod 
energy  to  short-period  energy  over  small  time  windows.  That  is  to  say,  the  more  complex  a 
body-wave  seismogram  the  higher  the  M f.  The  effect  on  the  Mjn^  diagram  would  be  that  tjv 
trend  line  for  complex  earthquakes  lies  to  the  high  Mg  side  of  the  overall  trend  line,  and  conse¬ 
quently  point  C  may  not  be  in  that  earthquake  population  with  a  comparable  complexity.  To  test 
tms  hypothesis,  complexity  as  a  function  of  Mg  and  mb  is  determined.  The  data  set  is  taken 
from  Lacoss.  Figure  IV-4  shows  the  Mg:mb  data  and  the  linear  trend,  Mg  =  1.6  mb  -  4.  M  m 
values  are  projected  onto  the  normal  to  the  linear  trend  and  plotted  against  Lacoss'  complexity  b 
values.  Thus,  comple.-*-  on  lines  parallel  to  the  overall  trend  line  is  determined.  The  result 
is  given  in  Fig.  IV-  5;  tlu  -shed  line  is  M;i  =  1.6  -  4.  The  data  show  Mm  pu  ifljrfk**.  R,f  i- 

groups  show  increasing  complexity  as  the  Mg:mb  trend  is  traversed  toward  increasing  values  of 


M3  along  the  perpendicular  to  the  overall  trend.  The  population  with  complexity  greater  than  5 
is  dominated  by  shallow  earthquakes  from  the  Kurile  Islands  leading  to  the  conclusion  that  the 
magnitude  of  the  phenomenon  is  a  function  of  source  location  and/or  mechanism.  For  either 
population,  the  result  is  that  the  trend  line  in  Fig.  IV- 3  moves  away  from  point  C  and  so  tends 
to  separate  the  Mg:mb  of  the  multiple  explosion  from  the  earthquake  population  with  comparable 
compi  y.  The  short-period  sum  trace  has  a  complexity  of  8,  and  the  corresponding  position 
of  po  t'  is  shown  in  Fig.  IV-5.  Thus,  the  total  effect  of  this  series  of  explosions  is  not 
earthquako-like  in  both  Mg:mb  and  complexity. 

To  simulate  a  pP  depth  phase  requires  producing  a  pulse  that  is  inverted  with  respect  to  an 
earlier  arrival.  While  it  is  impossible  to  do  it  exactly,  two  characteristics  of  single-explosion 
seismograms  written  on  typical  short-period  instruments  make  it  possible  to  produce  a  second¬ 
ary  pulse  whose  polarity  appears  inverted.  The  features  are  that  the  basic  pulse  is  essentially 
monochromatic  and  that  the  amplitude  of  tie  first  half  cycle  is  smaller  than  several  seconds  of 
the  early  coda.  Figure  IV-6  shows  five  of  the  same  source  functions  with  amplitude  ratios  of 
S.  3,  7.  4.  ->  anil  delays  of  0.  4,  4.7.  5.1.  7.8  sec  and  the  sum.  In  the  sum.  the  three  apparent 
arrivals  represent  P,  pP.  and  sP.  The  actual  first  motion  of  the  second  arrival  is  hidden  in 
the  coda  of  the  first.  Smaller  sources  are  advanced  by  one  cycle  and  delayed  by  one-half  cycle 
relative  to  the  second  major  arrival.  The  result  is  an  apparent  inversion.  The  quality  of  the 
•  aversion  compared  to  real  data  is  good.  The  Rayleigh  wave,  shown  in  Fig.  1V-7.  is  not  signif¬ 
icantly  enhanced  in  magnitude  and  so  the  M  -irv  value  for  the  sum  is  about  the  same  as  for  either 
oi  the  two  large  events.  Such  a  situation  does  not  necessarily  enable  discrimination,  as  problem 
events  with  similar  waveforms  exist.  The  success  of  the  above  simulation  critically  depends  on 
the  essentially  monochromatic  nature  of  the  source  waveform.  A  broadband  recording  of  such  a 
series  of  explosions  would  present  an  entirely  different  situation. 

In  conclusion,  the  possibility  of  confusing  a  series  of  explosions  with  an  earthquake  does 
exist.  A  hether  discrimination  can  be  made  by  redefining  A1  .mv  as  a  function  of  complexity,  by 
developing  numerical  schemes  that  are  capable  of  detecting  closely  spaced  functions  in  a  seismo¬ 
gram,  or  bv  recording  broadband,  is  a  subject  for  future  study. 

T.  E.  Landers 

B.  A  STUDY  OF  LASA  EXPLOSION  DETEC  TION  CAPABILITY 

WITHIN  THE  CODAS  OF  ASIAN  EARTHQUAKES 

A  simple  scheme,  employed  to  estimate  the  LASA  explosion  detection  capability  within  the 
codas  of  earthquakes  from  certain  Asian  regions,  consisted  of  the  following  steps.  The  phased 
sum  or  beam  was  formed  using  the  LASA  short-period  traces  of  an  explosion  and  an  Asian  earth¬ 
quake.  In  both  cases,  the  beams  were  aimed  at  the  explosion  site.  A  scaled  20-sec  segment  of 
the  explosion  beam  was  then  repeatedly  added  to  the  earthquake  coda.  A  narrow  bandpass  filte: • 
and  then  a  standard  detection  process  were  applied  to  the  beam  sum.  The  number  of  explosions 
missed  by  the  detection  process,  as  a  function  of  the  earthquake-explosion  magnitude  difference, 
gave  the  results  of  the  experiment. 

With  this  study  to  date,  four  pairs  of  explosion-earthquake  locations  have  been  considered: 

fane  I  Explosions  -  presumed  Soviet  test  site  (STS)  near  Semipalatinsk 

in  Kazakh  S.S.  R. 

Earthquakes  -  Kurile  Islands  region 
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Case  II 


Case  III 


Case  IV 


Explosions  -  STS 

Earthquakes  —  Kamchatka  Peninsula  region 

Explosions  -  hypothetical  test  site  in  the  Kurile  Islands 
Earthquakes  —  Kurile  Islands  region 

Explosions  -  hypothetical  test  site  near  the  Kamchatka  Peninsula 
Earthquakes  —  Kamchatka  Peninsula 

Experiments  related  to  Cases  I  and  II  are  straightforward  since  recordings  of  presumed 
STS  explosions  and  Kurile  Islands  and  Kamchatka  earthquakes  exist  in  our  data  library.  We  have 
no  recordings  of  explosions  in  the  Kamchatka  or  Kuriles  region,  but  an  attempt  to  synthesize 
Cases  III  and  IV  has  been  made  by  adding  an  STS  explosion  beam  trace,  aimed  at  STS,  to  a 
Kurile  earthquakes  beam  when  the  earthquake  beam  was  formed  on  a  site  in  the  Kurile  Islands. 

An  actual  Kurile  explosion  may  look  quite  different  than  an  explosion  in  central  Asia,  of  course, 
but  the  ruse  is  needed  to  perform  the  experiments.  The  reason  that  the  experiments  of  Cases  III 
and  IV  are  desirable  is  that  they  are.  in  a  sense,  extreme  eases.  We  expect  earthquake  coda 

attenuation  due  to  beamforming  to  be  a  minimum  when  the  array  beam  is  aimed  near  to  the  earth¬ 
quake  epicenter. 

The  detection  scheme  used  here  was  modeled  after  that  used  by  the  Seismic  Array  Analysis 
Center  ISAAC)  and  is  described  in  Ref.  3.  The  method  compares  the  ratio  of  a  short-term  aver¬ 
age  (STA)  rectified  signal  to  a  similar  long-term  average  (I.TA)  and  declares  a  detection  when 
this  ratio  exceeds  a  threshold  for  a  certain  number  of  comparisons.  No  search  was  made  to 
find  an  optimum  coda  detection  process;  a  standard  process  was  used  with  threshold  parameters 
which  gave  relatively  few  false  alarms. 

The  STA  was  computed  using  the  following  formula  from  Ref.  t: 

P-1 


STA(nAt)  =  P_1  £  I B  f  ( n  —  p )  At  ]  | 


p=0 

Here.  B(nAt)  represents  the  filtered  beam  sampled  at  nAt.  In  practice.  At  was  set  at  0.2  sec 

and  P  =  3;  the  STA  was  computed  every  0.6  sec.  The  I.TA  was  computed  using  the  recursive 
relation 


l.TA 


2'”  STAn  +(1  -2a)  LTA  4 
n  n- 1 


Values  of  1  and  4  are  used  for  ,j  and  a.  respectively.  Because  the  effective  length  of  this  LTA 
was  greater  than  the  20-sec  interval  at  which  explosions  were  added,  the  LTA  was  computed  using 
the  filtered  earthquake  beam  to  which  no  explosions  had  been  added.  This  means  that  the  pre¬ 
viously  added  explosions  did  not  contaminate  the  LTA  upon  which  the  detections  were  based.  The 
STA  of  this  beam  was  also  computed  in  order  to  check  for  false  alarms.  If  the  STA/LTA  ratio 
exceeded  10  dB  for  three  consecutive  computations,  i.e.,  1.8  sec,  a  detection  was  declared.  An 
example  of  the  method,  as  applied  to  Case  III,  is  given  in  Fig.  IV-8. 

The  presumed  explosion  used  in  this  study  occurred  on  8  September  1968  and  was  assigned 
a  body-wave  magnitude  of  S.S  by  the  U.S.  Coast  and  Geodetic  Survey  (USCGS).  This  event  was 
scaled  to  ten  other  magnitudes,  down  to  3.4,  in  order  to  repre.,  mt  a  wide  range  in  explosion 
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size.  The  locations  of  the  explosions,  taken  for  beamforming  purposes,  and  the  regions  from 
which  the  earthquakes  were  taken  are  summarized  below: 

Case  Explosions  Earthquakes 

I  50°N,  oO'E  40-50°N,  145-155”E 

II  50°N,  8l°E  50-60°N,  155-170°E 

III  45°N,  150’  E  40-50-N,  145-155'E 

IV  SS-N,  160°E  50-60°N,  155-170“E 

The  detection  criterion  was  applied  to  the  earthquake  trace  (in  order  to  check  for  false 
alarms)  and  to  the  11  earthquake-plus-explosion  traces  for  each  of  the  four  cases  using  16  earth¬ 
quakes  from  the  Kurile  Islands  region  and  18  earthquakes  from  the  Kamchatka  region.  The  num¬ 
ber  of  explosions  missed  (N)  was  counted  in  each  experiment.  If  N  was  less  than  the  number 
of  explosions  added  before  the  "arthquake,  due  to  high  noise  conditions,  that  trace  was  disre¬ 
garded.  Otherwise,  the  number  of  explosions  missed  was  recorded  and  plotted  against  the  mag¬ 
nitude  difference  between  the  earthquake  and  scaled  explosion.  These  results  are  given  in 
Figs.  IV-9  through  IV- 12. 

Figure  IV-9  shows  the  number  of  explosions  missed,  in  the  conditions  of  Case  I,  plotted 
against  the  earthquake-explosion  difference,  -  mx.  Figures  IV-10  through  IV-12  give  anal¬ 
ogous  displays  for  Cases  II  through  IV,  respectively.  The  straight  lines  drawn  on  these  figures 
represent  the  author's  estimate  of  the  maximum  number  of  explosions  missed  as  a  function  of 
magnitude  difference.  Although  parameterized  versions  of  these  lines  will  be  used  in  the  follow¬ 
ing  section  to  extend  these  results,  the  nature  of  these  lines  and  the  distribution  of  points 
merit  some  discussion  at  this  time.  Cases  I  and  II  (Figs.  IV-9  and  -10)  are  similar  in  that  the 
explosions  are  from  the  same  site  and  the  eartnquakes  are  from  regions  contiguous  to  each  other. 
In  spite  of  this,  the  scatter  of  points  is  greater  in  Figs.  IV-10  and  -12  than  in  Figs.  IV-9  and  -11; 
this  probably  reflects  the  wider  geographic  distribution  of  the  earthquakes  used  in  Cases  II  and  IV. 
It  is  also  of  interest  to  compare  the  trend  of  the  points  in  Cases  I  and  II  with  that  of  Cases  III 
and  IV.  The  former,  although  they  have  a  lower  zero  intercept,  exhibit  a  steeper  slope  than  the 
latter.  The  slope  of  the  "maximum  number  missed"  lines  is  largely  controlled  by  the  number 
of  missed  detections  at  high  magnitude  differences.  This  leads  to  the  conclusion  that  the  .:oda 
from  large  events  becomes  increasingly  incoherent  with  increasing  time  after  the  first  arrival. 
Alternatively  stated,  the  effectiveness  of  beamforming  on  location  A  suppress  the  coda  of  an 
event  from  location  B  decreases  with  time  after  the  arrival  of  the  initial  P-wave  from  location  B. 

The  application  of  the  detection  scheme  to  the  earthquake  trace  alone  gave  rise  to  nine  false 
alarms  in  the  total  of  68  beams  processed  (two  beams  for  each  event  in  each  region).  Of  these 
spurious  detections,  five  were  due  to  the  surface  reflection  phase  pP,  two  were  due  to  a  small 
interfering  event,  and  the  remaining  two  were  genuine  false  alarms. 

An  attempt  to  extend  these  results  in  order  to  estimate  the  maximum  total  time  that  an  explo¬ 
sion  of  given  magnitude  in  the  Kurile  Islands,  near  the  Kamchatka  Peninsula,  or  at  STS  would 
have  gone  undetected  at  LASA  due  to  the  long-term  seismicity  of  the  Kuriles  or  Kamchatka  region 

is  given  in  the  next  section.  ,  „ 

J.  R.  Filson 

C.  ON  ESTIMATING  THE  EXPLOSION  MASK  TIME  AT  LASA 

DUE  TO  THE  SEISMICITY  OF  CERTAIN  REGIONS 

In  this  section,  an  attempt  is  made  to  extend  the  results  of  Sec.  B  in  an  effort  to  estinia’e 
the  detection  capability  of  LASA  over  an  extended  period  of  time  for  each  of  the  four  cases 
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previously  discussed.  Let  the  functions  M,  E.  and  t  be  defined  as  follows: 

T,  s2>  The  number  of  earthquakes  of  magnitude  m_  that 

occurred  in  the  region  s_  over  the  extended  period 
of  time  T  c 

h’(niq  “  ^2*  An  est*mate  of  the  time  that  an  explosion  of  magn:tude 

H  mx  at  site  sj  would  go  undetected  at  LASA  in  the 

P-wave  coda  of  an  earthquake  of  magnitude  m  from 
rer  on  s2  1 

T<"V  "r’  3i>  s2^  An  estilnate  of  the  total  time  that  an  explosion  of  mag¬ 

nitude  mx  at  site  Sj  would  have  gone  undetected  due 
to  the  natural  seismicity  of  region  s_  over  the  time 
period  T.  i 

The  assertion  is  made  that  the  total  mask  time  r  may  be  represented  as 

r(mx,  T,  Sj.Sj)  =  V  M(mqll  s2)  E(mq  -  mx>  Sj.  s2)  (IV-1) 

where  the  summation  is  over  m  . 

q 

The  results  of  Sec.  R  give  us  an  opportunity  to  estimate  E  and  thus  t  using  published 
seismicity  lists.  The  straight  lines  on  Figs.  1V-9  through  IV-12  represent  the  maximum  number 
of  explosions  (N),  spaced  at  20-sec  intervals,  that  would  be  missed  at  LASA  for  an  earthquake- 
explosion  magnitude  difference  mq  -  m^.  These  lines  are  of  the  form 


N(m  - 

q 

-  mx.  Sj,  s2>  =  6  exp[f  (mq  -  m^)) 

(IV- 2) 

where  the  parameters 

6  and  £  depend  on 

s.  and  s_.  The  values  of  these 

■  parameter.-,  for  each 

of  the  four  cases  are  given  below: 

Case 

6 

£ 

I 

STS 

Kuriles 

1  9 

1.66 

II 

STS 

Kamchatka 

2.9 

1.64 

III 

Kuriles 

Kuriles 

4.5 

1.20 

rv 

Kamchatka 

Kamchatka 

5.5 

1.26 

Once  N  has  beei.  evaluated  for  a  given  case  and  a  given  m  -  mx.  it  may  be  used  to  estimate  E  by 

E(mq  -  mx,  sr  s2)  ^  At  •  N(mq  -  mx.  sr  s2)  .  (IV- 3 ) 

Here.  At  is  the  time  interval  at  which  the  explosions  were  added  (20  sec)  and  E  now  is  an  '-st.- 
mate  of  the  maximum  mask  time  due  to  an  event  magnitude  difference  m  -  m  .  Equations  (1V-2) 
and  (IV- 3)  make  possible  the  use  of  Eq.  (IV-1  )  in  order  to  estimate  t. 

There  are  certain  qualifications  that  must  be  made  before  the  results  of  such  estimates  are 
discussed.  The  time  period  T  was  taken  as  the  decade  1963  >.o  1072,  and  the  USCGS  and  NOS 
seismicity  lists  for  this  period  were  used  to  determine  the  function  M.  Although  these  lists  mav 
lack  sensitivity  below  a  certain  magnitude,  they  represent  the  most  stable  source  of  world-wide 
seismicity  information  over  an  extended  period  of  time.  Figures  IV-1  3  and  -14  show  the  cumula¬ 
tive  seismicity  distribution  for  the  Kurile  Islands  and  the  Kamchatka  Peninsula  for  the  period 
1963  to  <972.  The  Kurile  Islands  region  was  taken  as  40-50”N  and  145-155"E,  the  Kamchatka 
region  50-60”N  and  195-170” E.  These  curves  exhibit  the  familiar  flattening  feature  at  low  mag¬ 
nitudes  due  to  the  insensitivity  of  the  location  nei  vork.  Ti,js  insensitivity  is  probably  significant 
below  m^  =4.5  in  both  regions. 
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This  effect  gives  rise  to  a  specific  qualification  on  the  use  of  Eqs.  (IV-1 )  through  (IV  ?)  and 
(IV- 3).  Because  of  the  form  given  Eq.  (IV-1),  the  values  assigned  to  its  parameters,  and  the  de¬ 
tection  scheme  used,  an  explosion  may  be  hidden  by  an  earthquake  of  a  smaller  magnitude.  Thus, 
in  order  to  estimate  the  total  time  that  a  magnitude  4.5  explosion  would  be  hidden,  the  total  num¬ 
ber  of  magnitude  4.3  earthquakes  is  required.  Since  the  seismicity  function  M  is  certainly  in¬ 
complete  at  magnitude  4.3,  an  estimate  of  r  for  mx  =  4.5  must  certainly  be  low.  The  severity  of 
this  effect  increases  with  the  value  of  the  parameter  6.  In  the  use  of  Eq.  (IV-2),  all  fractional 
detections  were  rounded  down  to  the  nearest  integer  before  using  N  in  Eq.  (IV- 3).  This  is  a  fairly 
safe  procedure  because  all  of  the  points  of  Figs.  IV-9  through  IV- 12  lie  below  the  straight  lines. 

It  is  also  a  necessary  procedure  since  large  negative  values  of  m  -  mx,  even  though  they  might 
give  rise  to  values  of  0  <  E  «  1  if  N  was  allowed  to  float,  appreciably  contribute  to  t  for  M 
large  and  mq  small.  In  this  case,  the  explosion  would  be  the  much  larger  arrival  on  the  record 
a  there  would  be  no  question  ot  a  much  smaller  earthquake  hiding  it.  In  the  light  of  the  above 
discussion,  the  condition  that 

6  exp[«(nV  -  mx))<  1.0  (IV_4) 

when  estimating  r  for  a  given  mx>  where  m^  is  the  magnitude  above  which  M(m  ,  T.  s 2)  is  be¬ 
lieved  to  be  complete,  assures  the  validity  of  r  in  the  context  of  the  experiment? 

Figures  IV -15  through  IV -1R  show  the  results  of  the  method  for  Cases  I  through  IV,  respec¬ 
tively.  The  estimates  of  -  in  these  figures  are  maximum  estimates  in  the  sense  that  the  param¬ 
eters  of  Eq.  (IV-2)  were  chosen  for  a  line  that  lay  rt,ve  all  the  data  points.  Thus,  the  ordinate 
in  these  figures  represents  the  maximum  estimate  of  total  time  in  minutes  during  the  10-year 
period  that  an  explosion  of  magnitude  mx>  at  the  site  indicated,  would  have  gone  undetected  due 
to  the  seismicity  of  either  the  Kurile  Islands  or  the  Kamchatka  region.  The  arrows  on  these  fig¬ 
ures  indicat?  an  mx  below  which  the  estimates  are  not  strictly  valid.  The  placing  of  these  arrows 
is  based  on  the  condition  of  Eq.  (IV- 4)  where  rrC  r  4.5.  This  implies  that  the  seismicity  data  are 
considered  complete  down  to  mq  ,  4.5.  if  the  reader  chooses  to  pick  an  alternate  m'  ,  these  ar¬ 
rows  will  be  shifted  to  the  right  or  left  of  the  number  of  magnitude  units  that  the  new^n'  differs 
from  4.5.  The  straight  solid  lines  on  these  figures  represent  an  approximation  of  the  linear  ex¬ 
tension  of  the  points  to  the  right  of  the  arrows,  to  mx  values  to  the  left  of  the  arrows  where  the 
plotted  points  are  considered  to  be  low. 

As  an  example,  these  figures  allow  us  to  estimate  the  maximum  total  time  that  an  explosion 
of  magnitude  5.5  would  have  gone  undetected  at  LASA  during  the  10- year  period  (1963  to  1972) 
for  each  of  the  four  cases.  The  results  are  summarized  below: 

Case  Total  Mask  Time 

I  205  minutes  (3.4  hours) 

11  135  minutes  (2.3  hours) 

111  Approximately  2000  minutes  (33.4  hours) 

,v  Approximately  1000  minutes  (16.7  hours) 

Again,  these  estimates  are,  in  a  sense,  maximum  ones  and  they  refer  to  only  one  array  at  which 
a  standard  detection  scheme,  not  specialized  for  coda  detection,  was  applied. 

J.  R.  Filson 
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Fig.  IV-1.  Sum  of  eight  identical  short-period  waveforms  to  simulate 
a  complex  seismogram.  1 
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r  lg.  tv-5,  I  omplexity  on  lines  parallel 
to  linear  trend  (solid  line  in  Fig.  IV-4). 
Dashed  line  is  IV1  =  1.6  m,  -  4. 


Fig.  IV-9.  Points  represent  numbe 
aimed  at  STS,  due  to  Kurile  Island; 
nitude  difference  (mq  -  mx)  between 
is  relevant  to  Case  I  described 

maximum  number  of  explosi - 

sions  were  added  at  20-sec  intervals 


r  oi  .-jib  explosions  undetected  on  LASA  bean, 
i  earthquakes.  These  points  arc  plotted  vs  mag- 
earthquake  (niq )  and  explosion  (mx).  This  figure 
in  text.  Straight  line  represents  author's  estimate  of 
ions  missed  as  a  function  of  magnitude  difference  lixplo- 
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Fig.  IV-15.  Result  of  applying  Eq.  (IV-1) 
to  data  of  Figs.  IV-9  and  -10,  Case  I. 
Points  represent  estimate  of  maximum 
time  (in  minutes)  that  magnitude  mx  ex¬ 
plosion  at  STS  would  have  gone  undetected 
by  LASA  due  to  Kurile  Islands  seismicity 
during  decade  1963  to  1972.  Arrow  rep¬ 
resents  magnitude  below  which  results 
are  not  considered  strictly  valid.  Line 
represents  approximate  linear  extension 
of  points  at  magnitudes  greater  than  that 
of  arrow. 


Fig.  IV-16.  Diagram  analogous  to  Fig.  IV-15 
and  relevant  to  Case  II.  Points  represent  es¬ 
timate  of  time  (in  minutes)  that  magnitude  mx 
explosion  would  have  gone  undetected  by  LASA 
due  to  Kamchatka  seismicity  during  decade 
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Fig.  IV- 17.  Diagram  analogous  to  Fig.  IV- 15 
and  relevant  to  Case  III.  Points  represent 
estimate  of  time  (in  minutes)  that  hypotheti¬ 
cal,  magnitude  mx  explosion  in  Kurile  Islands 
would  have  gone  undetected  by  LASA  due  to 
Kurile  Islands  seismicity  during  decade  1963 
to  1972. 
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Fig.  IV-18.  Diagram  analogous  to  Fig.  IV-15 
and  relevant  to  Case  IV.  Points  represent 
estimate  of  time  (in  minutes) that  hypothetical 
explosion  of  magnitude  mx  near  Kamchatka 
Peninsula  would  have  gone  undetected  by 
LASA  due  to  Kamchatka  seismicity  during 
decade  1963  to  1972. 
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V.  GENERAL  SEISMOLOGY 


OBSERVATIONS  OF  PcP  AND  P  PHASES  AT  LASA  AT  DISTANCES 
FROM  40°  to  84° 

In  the  SATS  dated  31  December  1972  (DDC  AD-757360),  we  reported  observations  of  P  and 
PcP  amplitudes  in  the  distance  range  26°  to  40°  from  LASA.  Here,  we  extend  the  range  of  obser¬ 
vation  with  data  from  40°  to  84°  distance  with  depths  from  0  to  600  km.  Included  in  the  population 
are  the  explosion  LONGSHOT,  and  four  presumed  Soviet  explosions  -  one  in  the  Ural  mountains, 
two  from  Novaya  Zemlya,  and  one  from  Eastern  Kazakh. 

In  studies  of  the  core-mantle  boundary,  one  of  the  perplexing  features  of  short-period  data 
has  been  that  observed  PcP/P  ratios  are  strongly  scattered,  with  average  values  several  times 
larger  than  predicted  by  reasonable  earth  models.  There  are  many  possible  causes  for  this 
phenomenon,  summarized  by  Vinnik  and  Dashkov,1  but  one  popular  explanation  is  that  the  earth 
has  low  Q  values  for  P-waves.  Kan  .mori2  found  that  some  of  the  scatter  for  events  of  different 
depth  could  be  reduced  by  assuming  a  highly  attenuative  earth  with  Q  values  for  short-period 
compressional  waves  which  are  comparable  to  Q  values  for  longer-period  shear  waves.  Ibrahim3 
studied  P  and  PcP  data  from  explosions  only,  and  deduced  that  a  Q  model  similar  to  Kanamori's 
would  be  required  to  reduce  the  data  scatter. 

This  has  created  somewhat  of  a  paradox  because  short-period  attenuation  studies  by  Asada 
and  Takano,  Archanibeau  et  aL,  Frasier  and  Filson/*  and  others  show  that  Q  values  for  short- 
period  P-.vaves  are  an  order-of-magnitude  higher  than  those  deduced  by  Kanamori. 

We  present  LASA  data  here  showing  that  PcP/P  ratios  are  a  strong  function  of  focus  depth. 
However,  we  also  demonstrate  that  low  Q  models  do  not  adequately  resolve  the  problem  and  that 
the  scatter  in  PcP/P  ratios  is  due  to  amplitude  fluctuations  of  P  rather  than  PcP  waves. 

Figure  V-l  shows  the  distribution  of  data  in  depth  and  distance  from  LASA.  The  open  circles 
are  events  to  the  southeast  of  LASA  within  an  azimuthal  sector  from  131°  to  1  51  ° .  The  solid 
circles  are  events  to  the  northwest  within  the  sector  303°  to  319°.  Also  shown  by  x's  are  the 
explosion  I.ONGSHOT  and  four  presumed  Soviet  explosions.  The  Soviet  events  all  have  azimuths 
nearly  due  north  from  LASA. 

The  steered  beams  for  about  half  the  data  used  in  this  study  are  displayed  in  Fig.  Y-2.  Each 
pair  of  traces  consists  of  a  P  beam  above  the  corresponding  PcP  beam,  for  the  same  event, 
normalized  for  display  put  poses  to  approximately  the  saMe  gain.  he  events  are  arranged  ver¬ 
tically  in  order  of  increasing  distance.  Those  events  marked  S  are  from  the  southeast  sector 
from  LASA,  all  others  being  from  the  northwest.  By  steering  beams  for  each  phase,  the  receiver 
crustal  vari-tions  from  sensor  to  sensor  are  minimized  and  the  noise  is  reduced. 

In  tig.  \  -3,  the  PcP/P  ratios  of  A/T  for  the  LASA  beams  are  shown  vs  epicentral  distance. 
The  curve  is  the  theoretical  PcP/P  ratio  predicted  by  Johnson's  CIT  208  earth  velocity  model,7 
assuming  the  B1  core-mantle  model  of  Bullen  and  Haddon.^  No  earth  attenuation  has  been  in¬ 
cluded  in  the  model.  Because  of  this,  the  theoretical  curve  varies  insignificantly  with  depth  of 
focus,  and  may  be  taken  to  represent  all  depths. 

A  noticeable  feature  of  the  data  is  the  discontinuity  around  60°,  where  most  of  the  PcP/P 
data  suddenly  drop  from  above  to  below  the  curve.  This  is  clearly  a  depth  effect.  Almost  all 
the  data  above  the  curve  for  delta  less  than  60°  have  shallow  depths  less  than  200  km.  Data  below 
the  curve  with  delta  greater  than  60°  are  almost  all  deep  events.  The  five  explosions  are  all 
above  the  curve.  The  earthquake  at  68°,  nearly  coincident  with  the  explosion,  has  an  NOS  depth 
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fixed  at  33  km.  At  58°,  the  solid  circle  below  the  curve  has  a  depth  of  393  km,  and  has  a  PeP/P 
ratio  several  times  lower  than  the  shallow  events  at  the  same  delta. 

Ihis  effect  is  more  clearly  seen  in  Fig.  V-4  where  the  PcP/P  data  are  plotted  vs  depth  of 
event.  Here  we  see  that  the  explosion  data  at  zero  depth  are  grouped  with  shallow  earthquakes, 
with  PcP/P  ratios  3  to  4  times  those  for  the  deep  events.  Figure  V-5  shows  the  data  of  Fig.  V-4 
after  normalizing  all  ratios  to  a  distance  of  60°.  This  was  done  by  removing  the  trend  of  the 
theoretical  curve  in  Fig.  V-3  from  the  data  of  Fig.  V-4.  Thus,  PeP/P  ratios  at  distances  greater 
than  60°  were  amplified  to  remove  the  effect  of  spherieal  spreading  and  eore  reflection.  One 
sees  that  the  scatter  of  PcP/P  ratios  in  Fig.  V-5  is  less  than  in  Fig.  V-4,  and  that  this  depth  de¬ 
pendence  is  clearly  in  the  data. 

If  we  now  correct  the  observed  PcP/P  data  for  Kanamori's  Q  model,  that  is,  remove  the 
attenuation  effects,  we  obtain  the  data  shown  vs  distance  in  Fig,  V-6.  We  see  that  the  PcP/P 
ratios  have  been  lowered  everywhere  so  that  they  are  almost  all  below  the  theoretical  curve. 
However,  the  scatter  of  data  has  not  been  significantly  reduced.  The  PcP/P  ratios  at  40°  have 
been  severely  attenuated  and  the  explosions  are  still  significantly  above  the  deep  event  populations. 
These  PeP/P  ratios  are  normalized  to  a  delta  of  60°  and  plotted  vs  depth  in  Fig.  V-7;  this  figure 
shows  that  the  depth  effect  persists  in  spite  of  the  strong  attenuation  corrections  whieh  have  been 
made.  This  indicates  that  the  scatter  of  PeP/P  ratios  cannot  be  caused  by  a  laterally  homoge¬ 
neous,  low  Q  model  of  the  earth. 

In  Fig.  V -8,  the  PcP  amplitudes  are  shown  with  each  event  normalized  to  LASA  m(.  =  6.0, 
using  the  Gutenberg-Riehter  magnitude  formula.9  The  eurve  shows  the  heoretical  PcP  amplitude 
for  the  unattenuated  earth  model.  It  is  adjusted  in  amplitude  to  pass  through  the  data.  The  PcP 
amplitudes  show  no  separation  between  deep  and  shallow  events.  In  fact,  the  trend  of  data  from 
50°  to  83°  is  well-defined  with  less  scatter  than  seen  in  the  PcP/P  ratios. 

This  suggests  that  the  Gutenberg-Riehter  13  faetor,  which  strongly  depends  on  focal  depth 
for  all  distances,  applies  to  P  data  but  not  to  PeP  amplitudes.  That  is,  magnitude -6  events 
located  at  various  depths  produce  a  relatively  smooth  PcP  curve  without  any  further  depth 
corrections. 

Figure  V-9  shows  the  PeP  amplitudes  of  Fig.  V-8  plotted  as  a  function  of  depth  and  normalized 
to  A  =  60°  by  removing  the  trend  of  the  curve  in  Fig.  V-8.  We  see  that  the  clear  depth  trend  of 
the  data  has  disappeared,  and  the  explosion  data  at  zero  depth  are  more  consistent  with  the  deep 
focus  data. 

The  conclusion  to  be  drawn  from  these  data  is  that  the  scatter  of  PcP/P  ratios  at  LASA  is 
primarily  caused  by  the  depth  effeet  on  P-waves,  which  is  given  by  the  empirical  13  factor  of 
Gutenberg  and  Richter.  This  13  factor  cannot  be  produced  by  attenuation  and  spreading  of  rays 
in  laterally  homogeneous  earth  models.  Other  scattering  mechanisms  are  being  investigated  to 
explain  this  phenomenon. 

W.  Frasier 
IX  K.  Ohowdhury* 


*  Indiana  University. 
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Fig.  V-i.  Distribution  of  data  in  depth  and  distance  from  LASA.  Open  circles 
are  events  southeast  of  LASA;  solid  circles  are  events  northwest  of  LASA. 
Presumed  explosions  are  denoted  by  x's  and  are  all  northwest  to  north  of  LASA. 
These  symbols  are  used  on  all  subsequent  data  figures. 


Fig.  V-2.  P  and  PcP  beams  of  approximately  one-half  of  data  studied.  Invents 
marked  S  are  southeast  of  LA.SA,  all  others  being  northwest  of  LASA. 
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Am;m»i„aPCP/P/amplitude  rati°3  as  a  functio"  of  distance  from  l.ASA 
Amplitudes  are  n/sec  corrected  Tor  instrument  response.  Curve  is  theo¬ 
retical  ratio  expected  for  nonattenuative  earth  model  described  in  text. 
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Fig.  V-  4.  Pc  P/P  ratios  plotted  as  a  funct 


ion  of  depth  of  focus. 
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Fig.  V-7.  PcP/P  ratios  corrected  for  Kanamori's  Q  model 
and  normalized  to  A  -  60°. 
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Fig.  V-8.  PeP  amplitudes  (n/sec )  normalized  to  I,ASA  mb  =6.0  as  measured 
on  P-wavcs  for  each  event.  Curve  is  theoretical  PcP  amplitude  assuming  un- 
attenuated  earth  model. 
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Fig.  V-9.  PcP  amplitudes  of  Fig.  V-8  plotted  as  a  function  of  depth 
after  removing  trend  curve. 


